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rH ■ A numerical study is made of the spectra of a tight-binding hamiltonian on 

^ ■ square approximants of the quasiperiodic octagonal tiling. Tilings may be pure or 

random, with different degrees of phason disorder considered. The level statistics for 
k>( I the randomized tilings follow the predictions of random matrix theory, while for the 

perfect tilings a new type of level statistics is found. In this case, the first-, second- 
level spacing distributions are well described by lognormal laws with power law tails 
for large spacing. In addition, level spacing properties being related to properties of 
the density of states, the latter quantity is studied and the multifractal character of 
the spectral measure is exhibited. 
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Since the discovery of quasicrystalline alloys, many geometrical models have 
been built that are consistent with observed diffraction patterns. Tilings may have 
perfect quasiperiodic long range order with local symmetries that are not allowed for 
ordinary periodic structures, and which give rise to the five-fold, eight-fold or ten-fold 
symmetries of the diffraction patterns. Alternatively, one may consider randomized 
versions of perfect tilings. The approximants can be geometrically disordered to 
varying degrees by performing "phason" flips whereby a vertex is displaced and 
the tiling redrawn, resulting in a rearrangement of the tiles that surround that 
vertex, "allowed" permutations of the tiles. A typical random tiling (i.e. excluding 
some special cases that could be generated by the randomization procedure) has 
the same local symmetry as its perfect counterpart. The perfect tiling of course 
has exact symmetries such as its self-similarity property under inflation or deflation 
operations, while random tilings do not have this property. Their orientational order 
and the constraint imposed by the flxed shape of the tiles of the random tiling are 
sufficient, however, to give rise to sharp peaks in its diffraction pattern, similar to 
the Bragg peaks of the perfect tiling [|l| . Many existing phenomenological models to 
explain quasicrystal electronic properties are based on the construction of a pseudo 
Brillouin zone taking into account the brightest spots of the diffraction pattern, and 
at this level they would thus not distinguish between the perfect and random tiling 
models. The effect of phason disorder on electronic properties can however be very 
profound, as will be shown below. 

In this paper we make a numerical study of statistical properties of the eigen- 
value spectrum of a tight-binding Hamiltonian describing the hopping of electrons 
on vertices of a two dimensional tiling. Hopping is allowed between sites that are 
separated by one unit length (the tile side is taken to be of length one). The tilings 
are chosen to be square periodic approximants of the quasiperiodic octagonal tiling, 
containing upto 8119 sites (the largest approximant considered here). The hopping 
parameter is taken to be constant and set equal to unity, so that the resulting spec- 
tra are a pure consequence of the geometry, or site connectivity, of the tilings. The 
Hamiltonian is diagonalized numerically for each case, and the energy levels are used 
to compute the distribution P{s) of spacings of n^^ nearest-neighbor energy levels 
s = Ei^n — Ei (supposing the levels to be ordered). To place the current study in 
perspective, we mention previous work on electronic levels of a tight-binding model, 
namely studies of disordered periodic lattices on both sides of the metal-insulator 
transition 0. These show that in the metallic regime, the flrst-neighbor level spacing 
distributions P(s) have the forms of the Wigner distributions, obtained for gaussian 



random matrices |]^. In other words the energy levels of a real symmetric Hamil- 
tonian with random distributed onsite energies have the statistics of levels of the 
gaussian orthogonal ensemble (GOE), where P{s) ~ s^exp(— c/js^) where c^ is a 
known constant and (3 = 1. In the case of a magnetic flux traversing the lattice, the 
Hamiltonian belongs in the unitary class and levels obey gaussian unitary ensem- 
ble (GUE) statistics, with (3 = 2. Finally, in the insulating regime (strong disorder) 
P(s) has a Poissonian form for large spacings, P{s) ~ exp(— s). A second quantity of 
interest calculated for these crystalline models is the spectral rigidity T?{E)., which 
measures the fluctuation of the number of levels in a energy window of width E. 
T?{E) grows logarithmically with E at small E in the metallic diffusive regime 0, 
just as for the random matrix ensembles. This behavior is a consequence of corre- 
lations between levels, and gives a smaller T?{E) compared to that of uncorrelated 
levels (as occurs in the insulating regime). 

In our present work, we find that the randomized quasiperiodic tiling has level 
statistics of the GOE or GUE type, depending on the applied flux. The statistics of 
energy levels are thus similar to those of disordered tight binding models on crystals 
in the diffusive regime, even though in our case the disorder is purely geometric [Q. 

For the perfect approximants, we find a new type of level statistics. The perfect 
octagonal approximant has been studied previously by Benza and Sire [^ who pre- 
sented results on level statistics. They did not propose a fit to any functional form, 
noting principally that level repulsion occurs in the perfect quasi crystal, since P(s) 
drops to zero at small values of s. This is certainly true of the distribution that we 
propose in this paper. For the perfect case, the calculations show that the density 
of states (DOS) has huge fluctuations at all energy scales, indicating multifractal 
properties of the spectral measure. We have calculated the effective /(a) function 
of local singularity exponents a. Its maximum occurs aX a = Di = f{Di) (the most 
probable dimension) and has the maximum value fmax = Dp ~ 1 (the fractal di- 
mension). The last is an indication that the spectrum is gapless, in accord with the 
finding of 0. The calculated first-neighbor distribution P(s), as well as the second- 
and third-neighbor spacing distributions P^'^\s) and P^^\s) are well described by 
lognormal laws with power law tails at large spacings. To complicate matters, how- 
ever, the low moments of the spacing distributions apparently do not have the size 
dependence that would be indicated by the calculated /(a), while the high moments 
have a size dependence that is characteristic of a power law behavior s'"^ in the large 
s tail of P(s). It should be noted that the systems studied are comparatively small, 
and studying the next size of approximant will help in determining better the size 



dependence of these distributions. 

We now describe more precisely our calculations. The perfect periodic approx- 
imant is generated using the method described by Duneau et al 0. The largest of 
the three sizes of tiling used in these computations is the 5/4 approximant which 
contains 8119 atoms in the unit cell. The random phason disordered approximants 
are generated by an iterative process of displacing sites and computing the new links 
between sites in each cycle. The Hamiltonian for a given tiling is defined by the 
following site-projected form 

where the sum is taken over sites j linked to site i. The number of links emanating 
from a given site i varies between 3 and 8 in the perfect case, corresponding to the 
six local environments |0 of the octagonal tiling. In the randomized case new local 
configurations appear, although the mean coordination is z = 4 in both tilings. To 
reiterate, in this model the quasiperiodicity as well as the disorder are present purely 
in the connectivity and are of purely geometric nature. This may be compared with 
models in which one takes an explicit parameter in the Hamiltonian to generate the 
quasiperiodicity or disorder, such as a site-dependent potential. This distinction 
is not possible in one dimension, where the quasi-periodic nature of connections 
between sites must be coded by introducing for example two kinds of bonds, as in the 
Fibonacci chain |^. Periodic continuation of the LxL square tilings allows to impose 
boundary conditions of the form ipi^x + L,y) = e^'^ipi^x^y)] ^J^x^y + L) = ip{x,y). 
This is equivalent to adding a magnetic flux (p along the y or z axis. It is thus 
possible to eventually induce a GOE to GUE transition by varying (p. 

We discuss first the results for the geometrically disordered case. Their anal- 
ysis and interpretation is more straightforward and so will serve as a guide to the 
case of the perfect approximant. The numerical results obtained for P{s) (taking 
a normalized spacing s = (-Ej+i — Ei)/{W/N) where W is the band width and N 
the number of levels) and for ^"^{E) are shown in Figs, la and lb. The points 
are obtained for flux values = and = 1/2 , while the curves give the fit to 
the corresponding analytical expressions from random matrix theory (RMT). It has 
been shown, for weakly disordered metals that the RMT law for the spectral stiff- 
ness can be derived, both in perturbation theory f^ , and by means of a semiclassical 
model [|10] by assuming a diffusive motion for the particle. Both derivations obtain 
that the spectral stiffness follows the RMT logarithmic law in the low energy region 
T.^'^^E) oc -^\n(E) + est, crossing over to a linear behavior (in two dimensions). 



Our calculations yield a similar behavior of the spectral rigidity. This may indicate 
that the disordered quasicrystal possesses a diffusive regime, as for the disordered 
crystalline metals. More precisely, if the wavefunctions are localized in the quasiperi- 
odic tiling, that the localisation length must be greater than the system size, since 
on the scale of our samples the results are in good accord with the Wigner distri- 
butions. We note that we have not been able to disorder (in the geometric sense) 
the tilings sufficiently to see a crossover to strong localization. Such a crossover can 
be seen in tilings that have a strong energetic disorder (see note in |Q) wherein we 
take a onsite energy randomly from an energy interval of width comparable to the 
bandwidth of the perfect tiling. In the strongly localized regime, P{s) tends to the 
exponentially decaying Poissonian form 0. Our results indicate that in all likeli- 
hood, the geometrically randomized tilings still possess sufficient long range order 
that wavefunctions remain "extended" in the sense of not decaying exponentially, 
however, this remains to be established. 

In contrast to the random case where the fluctuations of the DOS are small 
compared to Poissonian, the perfect approximant, like one dimensional quasiperiodic 
systems such as the Fibonacci chain, has a self-similar DOS with huge fluctuations 
at all energy scales |jll|]. Quantities which characterize these multifractal properties 



are the generalized fractal dimensions D{q). To flnd these we consider a partition 
of the bandwidth W into M distinct boxes Si of probability Pi and size li where 
Pi = j^ and Hi is the number of level in the box of size li. The dimensions D{q) 
satisfy the following condition ([|T^) 

M 9-1 
1=1 'i 

with maxi(li) -^ and r(g) = {q — l)-D(g). In the simplest case, one takes a par- 
tition of boxes with equal size / = W/M. In that case, the measure pi is assumed 
to scale as Pi ~ /"* when / tends to zero and the number of aj between a and 
a + da is assumed to vary as Ni{a) ~ l~^^"^da. The function r(g) vs q and the 
effective /(a) vs a are related by means of a Legendre transform. Conversely, one 
can flx the weights of the boxes to be identical, allowing their widths to fluctuate, 
and calculate q^r). This is tantamount to calculating moments of spacings with 
respect to the underlying distributions P^")(s). The n*'' neighbor level spacing mo- 
ments li = s^ = {Ei^n — Ei)/W,i = 1,N are obtained by taking the partition of 
flxed probability p = n/N. 

We have computed the /(a) function using the two methods described above. 



In agreement with the earher work of Benza and Sire who found a single-band spec- 
trum for the pure hopping model, we find the maximum of /(a) to have the value 
Dp ~ 1 indicating the spectrum is gapless. Another numerical argument in favor of 
a spectrum with no (finite measure) gap is our finding that the biggest spacing Smax 
scales as ~ N^^'"^ {smax -^ for increasing A^ is a necessary condition to apply the 
second of the two methods described above). The information dimension defined 
earlier is Di = 0.98 and a fit to a parabolic form of /(a) close to its maximum 
gives its curvature l//i, /i = 0.03. Now a parabolic shape for /(a) can be shown to 
lead to lognormal spacing distributions for each of the nth nearest neighbor spacings 
s" (for n << N). The parameters of these distributions are related to Dj and /i. 
The distributions depend, as well, upon the system size. We have found, however, 
that when the distribution P{s) of rescaled variables s = N{Ei^i — Ei/W) is plot- 
ted, a single curve is obtained for the three sizes studied: A^ =239,1393 and 8119 
sites. The same remark applies to the two other distributions calculated, p(^'^). A 
comment on the calculation of P{s): it is well known that before computing any 
statistical spectral properties, the spectrum has to be unfolded f^. This procedure 
can be understood within a semiclassical scheme as being equivalent to taking only 
the fluctuating part of the integrated density of states. One keeps only the quan- 
tum interference corrections around the Thomas-Fermi zeroth h order term of the 
integrated density of states. However to have an efficient unfolding, the fluctuating 
part has to be a small perturbation around the zeroth h order term and the zero 
order term itself should not be strongly fluctuating. In weakly disordered crystalline 
systems this is the case. In contrast, the density of state of a quasiperiodic system 
like ours has huge fluctuations at all energy scales. One may consider that in this 
case at least one of the hypothesis breaks. Thus the unfolding procedure is not well 
deflned for our spectrum. We have therefore computed P{s) without any unfolding. 
The distribution of the logarithm of s is well fltted by a gaussian, so that P{s) 
is a lognormal function of the form: 

/ N 1 (ln(s) — In(sn))'^ , , 

ns) = -^^exp - ^ ^ '^^'" (3) 

with ln{so) = —B/2 for variables normalized to unity so that J^ sP{s)ds = 1. 
Level repulsion clearly occurs, as P{s) tends to zero with diminishing s at small s. 
The lognormal form is obeyed for spacings around the most probable value. At small 
s the repulsion seems linear rather than exponential as given by the lognormal, 
while at large spacings, s > 3W/N, the error function J^ P{s)ds behaves as ~ s~^ 
implying P{s) ~ s~^ (this agrees with the observed variation in the largest spacing 



as a function of the system size and with P{s) ~ 5^(1+"™"^)). A final observation on 
the dependence of P{s) on the flux traversing the tihng: we have done a separate 
calculation for each value of the flux 0. In the perfect approximants we find P{s) 
to be independent of the flux, unlike the case of weakly localized systems |^. 

In conclusion, we have numerical evidence for fundamental differences under- 
lying spectra of perfect and random quasicrystal approximants, the latter resem- 
bling disordered crystalline systems in so far as their level statistics are concerned. 
The perfect approximants have lognormal distributions for the level spacings, with 
crossover to a power law in the tails. Results for different sample sizes can be made 
to coincide by plotting distributions of suitably scaled variables. This result seems 
inconsistent with the size- dependence expected of a straight-forward multifractal 
energy spectrum and calculations on bigger systems should help resolve the issue. 
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Figure Captions 

la. Level spacing distribution for the randomized approximants; in the GOE case 

(0 = circles) and GUE (0 = 1/2 stars). Continuous lines show the corresponding 

Wigner distributions. 

lb. Spectral rigidity in the same two cases along with the RMT fit ( filled circles: 

GOE, open circles: GUE). 

2. f{a) distribution of singularity exponents a along with a fit to a parabolic shape. 

3a. Level spacing distribution for the perfect approximant fit to a lognormal form. 

In the inset the distribution of the logarithm of s is fit to a gaussian. 

3b. and 3c. Same as 3a for the next nearest and third neighbor spacing distribution. 



